Turbulence-condensate interaction in two dimensions 
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We present experimental results on turbulence generated in thin fluid layers in the presence of 
a large-scale coherent flow, or a spectral condensate. It is shown that the condensate modifies the 
third-order velocity moment in a much wider interval of scales than the second one. The modification 
may include the change of sign of the third moment in the inverse cascade. This observation may 
help resolve a controversy on the energy flux in mesoscale atmospheric turbulence (10-500 km): to 
recover a correct energy flux from the third velocity moment one needs first to subtract the coherent 
flow. We find that the condensate also increases the velocity flatness. 
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More often than not, turbulence coexists with a flow 
coherent across the system size. Such flows can be ex- 
ternally generated or appear as a result of spectral con- 
densation due to an inverse turbulent cascade (see e.g. 
[TJ HJ [31 21 [5J [B]). Understanding interaction between 
turbulence and a mean flow is of prime importance for 
many problems in astrophysics, geophysics, plasma con- 
finement, etc. From a fundamental viewpoint, it is inter- 
esting to understand how spectral condensation breaks 
symmetries that emerge in inverse cascades [7], and how 
the condensate suppresses turbulence level and turbulent 
fluctuations [SJ [51 H5]. From a practical viewpoint, at- 
mospheric physics presents arguably the most important 
cases of turbulence which interacts with the large-scale 
flows. 

Atmospheric motions are powered by gradients of so- 
lar heating. Vertical gradients cause thermal convec- 
tion on the scale of the troposphere depth (less than 
10 km). Horizontal gradients excite motions on the 
scales from 1000 to 10000 km. Both inputs are redis- 
tributed over wide spectral intervals by nonlinear inter- 
actions [lOl [TTJ [12] . The wavenumber spectra measured 
in the upper troposphere and in the lower stratosphere 
have shown two power laws: E(k) oc fc~ 5 / 3 for the scales 
between 10 and 500 km, and a steeper spectrum with 
E(k) cx fc~ 3 in the range (500-3000) km (similar to the 
spectra in Fig. 2) jT3j. Interestingly, such spectra appear 
in the Kraichnan theory of 2D turbulence [10], where 
E(k) = Ce 2 / 3 fc -5 / 3 corresponds to an inverse energy cas- 
cade and E(k) = C q r] 2 ^ 3 k~ 3 to a direct vorticity cascade, 
e and rj being the dissipation rates of energy and enstro- 
phy respectively. This prompted a two-source picture 
of atmospheric turbulence with a planetary-scale source 
of vorticity and depth-scale source of energy, where the 
large-scale spectrum is due to a direct vorticity cascade 
|14j . Alternatively, that spectrum can result from an in- 
verse cascade of inertio-gravity waves [3]. Yet another 
possibility is that the main energy at large scales is actu- 
ally not in turbulent fluctuations but in a long-correlated 
flow (condensate), which can be either generated by ex- 



ternal forces or appear in the process of spectral con- 
densation (turbulent counterpart of Bose-Einstein con- 
densation), as suggested in [4] and demonstrated exper- 
imentally here. No less controversial is the nature of 
the mesoscale 5/3-spectrum. Is it an energy cascade and 
what is the flux direction? 

In homogeneous turbulence, spectral energy flux is ex- 
pressed via the third-order moment of the velocity |15j : 
e = S 3 /r, where S 3 = [(((SV L ) 3 ) + (SV L (SV T ) 2 ))] /2. 
Here SVl and SVr are respectively longitudinal (l) and 
transverse (?) components of the velocity difference be- 
tween points separated by r. Angular brackets denote 
time averaging. Positive S3 corresponds to the inverse en- 
ergy cascade (from small to large scales). Measurements 
of 53 in the atmosphere gave a negative value in the in- 
terval 10-100 km, which was interpreted as the signature 
of the forward energy cascade [T5] . Here we demonstrate 
experimentally that a negative small-scale S3 in a system 
with an inverse cascade can be caused by a large-scale 
shear flow. 

Let us first consider how small- and large-scale parts 
of the velocity difference (respectively Sv and SV) con- 
tribute to the second and third velocity moments. The 
large-scale flow is spatially smooth so that SV ~ sr 
where s = V/L s is a large-scale velocity gradient and 
L s is the shear scale which depends on the system size 
and on the topology of the flow. Comparing ((SV) 2 ) = 
s 2 r 2 with ((Sv) 2 ) = C(er) 2 / 3 we see that the small- 
scale (turbulent) part dominates at the scales smaller 
than l t = C 3 / 4 s _3 / 2 e 1 / 2 . For the third moment, we 
compare ((Sv) 3 ) = er with the cross-correlation term 
(SV(Sv) 2 ) = srC(er) 2 / 3 and observe that the influence of 
SV extends to a much smaller scale Z„ = C~ 3 / 2 s~ 3 l 2 e 1 / 2 , 
because the dimensionless constant C is substantially 
larger than unity, as discussed below. At the scale Z*, 
((SV) 3 ) ~ (sr) 3 <C er. That attests to a non-locality 
of condensate-turbulence interaction: not only conden- 
sate breaks scale-invariance of turbulence, but it imposes 
different scales on different moments. 

The above estimates are true for a large-scale part pro- 
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duced by any source. In particular, when it is produced 
by an inverse cascade (as in the experiments described 
below) one estimates s as follows. Let the linear damp- 
ing rate a be smaller than the inverse turn-over time 
£<i/2 £ i/3£-3/2 £ or vortices comparable to the system 
size L. Then the flow coherent over the system size (the 
condensate) appears [2J [5J El [TTJ [HI [TJ5] with the veloc- 
ity estimated from the energy balance, aV 2 ss 2e, which 
gives s = V/L s = L~ X yj2e/a and 

fc t = 7r// t -^ s - 3 / 2 (Ca/2)- 3 /V/ 4 . (1) 

Note that this is not the condition that the turnover time 
at /t is a -1 , as in [T7]; incidentally Eq. [I] gives a correct 
estimate (kt = 1) for their conditions. The spectrum 
E(k) cx fc -3 at k < k t is due to the condensate [51 151 ITT], 
while E(k) = Ce 2/3 fc~ 5 / 3 is expected at k t < k < kf. 

Here we report the experiment in which the strength 
and the spectral extent of the condensate are varied by 
changing either a or L. The experimental setup shown in 
Figure[I]is similar to those described in [6j[T9j[20] but has 
a substantially larger number of forcing vortices (up to 
900), higher spatial resolution and larger scale separation 
(L/lf ss 30). Turbulence is generated electromagneti- 
cally in stratified thin fluid layers whose thicknesses are 
varied to achieve different a. A heavier non-conducting 
fluid (Fluorinert, specific gravity SG— 1.8) is placed at 
the bottom. A lighter conducting fluid, NaCl water so- 
lution (SG= 1.03), is placed on top. Non-uniform mag- 
netic field is produced by a square matrix of 30 x 30 
permanent magnets (10 mm apart). The electric cur- 
rent flowing through the top (conducting) layer produces 
J x B-driven vortices which generate turbulence. Square 
boundaries with L = (0.09 — 0.24) m are used. To visu- 
alize the flow, imaging particles are suspended in the top 
layer and are illuminated by a 1 mm laser sheet parallel to 
the fluid surface. Laser light scattered by the particles is 
filmed from above using a 12 Mpixel camera. Green and 
blue lasers are pulsed for 20 ms consecutively with a delay 
of (20-150) ms. In each camera frame, two laser pulses 
produce a pair of images (green and blue) for each parti- 
cle. The frame images are then split into a pair of images 
according to the colour. The velocity fields are obtained 
from these pairs of images using particle image velocime- 
try. The damping rate (in the range a — 0.05 — 0.5 s _1 ) 
is estimated from the decay of the total kinetic energy 
after switching off the forcing: E cx e~ at . 

Figure [2] shows the energy spectrum measured for large 
L = 0.235 m and an intermediate a — 0.16 s . The 
force wavenumber is kf ss 400m~ 1 ; E(k) cx k~ 3 at 
k > kf , while E(k) cx fc~ 5 / 3 at k t < k < kf. At 
k < k t ~ 80 m _1 , in the condensate range, the spectrum 
is steeper and close to k~ 3 . Due to the condensate, the 
spectrum has k~ 3 and fc -5 / 3 ranges for the large and in- 
termediate scales respectively, similarly to the Nastrom- 
Gage spectrum [13] . Spectra for different L and a are 
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FIG. 1: Experimental setup 
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FIG. 2: Kinetic energy spectrum measured for the largest 
box L = 0.235 m and intermediate damping a = 0.16s -1 . 
The guide lines show the power laws for different ranges: fc~ 3 
vorticity cascade, k~ 5 ^ 3 energy cascade and fc~ 3 condensate. 

shown in Fig. [3] At fixed a — 0.15 s _1 , the knee of the 
spectrum shifts from k t ~ 80 m _1 for L = 0.235 m to 
k t m 135 m _1 for L = 0.15m, Fig. [3^a). For fixed L, 
linear damping affects k t as shown in Fig. gb). Going 
from a = 0.15 s _1 to a = 0.06 s^ 1 , changes kt from 
80 m _1 to kt ~ 130 m -1 . These observations are in 
a good qualitative agreement with (JT|). By further re- 
ducing L, we achieve a regime when k t ~ kf, and the 
ran g e disappears, such that the entire spectrum is 
E(k) cx fc -3 , both above and below kf, as in There- 
fore, we can control the shape of the spectrum and the 
relative strength of the condensate with respect to tur- 
bulence. 

We now analyze two regimes, a weak and a strong con- 
densates, whose spectra are shown in Fig.[4][a,c). A weak 
condensate of Fig. [ija) was generated at a = 0.3 s _1 and 
L = 0.235 m, while a stronger condensate of Fig.[ijc) was 
obtained at a = 0.15 s _1 and L = 0.15 m. For the weak 
condensate case, Fig.BFb) shows normalized velocity mo- 
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ments. Skewness, Sk — Ssl/S 2L (open triangles) is pos- 
itive and it is in the range of Sk = 0-^0.2 (skewness of the 
Gaussian process is Sk = 0). Both longitudinal S^l = 
((SVl) 3 } and transverse S3T = (SVl(SVt) 2 } moments are 
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FIG. 3: Kinetic energy spectra (a) for different box sizes at 
a — 0.15 s , and (b) for different linear damping rates at 
L = 0.235 m. 

positive in the entire range of scales, in agreement with 
expectations for the inverse energy cascade. The Sk(r) 
dependence has a knee at about r — 0.04 m, which cor- 
responds to the knee in the spectrum at k t as 80 m -1 , 
Fig. ga). The Kolmo gorov constant is determined as 
C = E(k)e~ 2 / 3 k 5 / 3 , where e = S$/r. In the weak con- 
densate, at r < n/kt ~ 0.04 m, we have C = 5.6, which 
is close to the values C = 5.8 4- 7 previously obtained 
in numerical simulations of 2D turbulence (see [55] and 
references therein). At larger separations, r > ir/k t , the 
function Sk(r) grows fast and cannot be approximated 
by a constant. 

The third-order velocity moment differs markedly in 
the presence of strong condensate. Fig. |4jd) shows that 
skewness (open triangles) varies in this case in a much 
wider range, Sk = —0.2 4 0.5. Flatness, F — S±lIS\ l 
(open squares), also varies, in the range F — 2.545, while 
its Gaussian value is F = 3. For the stronger condensate, 
the spectrum scales as E(k) cx k~ 5 / 3 in the range kf > 
k > k t ~ 125 m . S%(r) and Sk change sign at r t = tt jk t 
(Fig.EFd)). Such an S3 dependence on r resembles ^(r) 
measured in the lower stratosphere [16]. Note that in our 
case all the driving comes from small scales and there 
is no direct cascade at all, yet S3 is strongly modified 
compared with the weak condensate case. 

The generation of the mean flow can be revealed by a 
temporal averaging of the instantaneous velocity fields: 
V(x,y) = (1/iV) Yl n =i V( x ' Vi The power spectrum 
of the mean flow [5J [5] is close to E(k) cx k~ 3 . The ve- 
locity field contains both the mean component and tur- 
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FIG. 4: (a,b) Weak condensate case, (c-f) Strong condensate, 
(a) and (c) are the kinetic energy spectra, (b) Normalized ve- 
locity moments for weak condensate: skewness Sk (triangles) 
and flatness F (squares). The moments are computed be- 
fore (open symbols) and after (solid symbols) the mean flow 
is subtracted from individual velocity fields, (d) Skewness 
(triangles) and flatness (squares) for the case of strong con- 
densate before subtracting the mean flow. After subtracting 
mean flow in the stronger condensate case: kinetic energy 
spectrum (e), (f) skewness (triangles) and flatness (squares). 

bulent velocity fluctuations: SV = SV + SV. From the 
data corresponding to Fig. |4|c) (strong condensate) , we 
estimate that ((SV) 2 ) differs from ((SV) 2 ) by about 20- 
30%. However ((SV) 3 ) and ((SV) 3 ) differ by orders of 
magnitude and even the sign can be different. Note that 
signs, values and functional dependencies S^(r) vary a lot 
for different topologies of the condensate flows and also 
depend on the mean shear in such a flow. 

To recover the statistical moments of the turbulent ve- 
locity fluctuations we take N = 350 instantaneous veloc- 
ity fields, subtract their mean flow and then compute the 
Fourier spectrum and the structure functions. The result 
for the weak condensate case is shown as solid symbols 
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in Fig. |4|b). Mean subtraction brings skewness (solid 
triangles) much closer to its Gaussian value of S = 0, 
while flatness (solid squares) is only slightly higher than 
its value in isotropic turbulence, F = 3 [T5] . 

The result for the strong condensate is shown in 
Figs. |4je,f ) . The subtraction of the mean restores the 
^-5/3 ran gg p j n f ac t ; the E(k) scatter in Fig.Qe) is less 
than in the total spectrum of Fig. |4jc) . The subtraction 
has even more dramatic effect on the higher moments. As 
seen in Fig. |4jf), there is less variability in both skewness 
and flatness. S3 is now positive and it is a linear function 
of r in the "turbulence" range, at r < r t . The spectral 
energy flux is deduced as e = S^/r. The value of the 
Kolmogorov constant C = E(k)e~ 2 / 3 k 5 / 3 w 7.7 appears 
to be slightly higher than in the weak condensate case, 
but is still close to the values obtained in numerical sim- 
ulations. The recovery of the linear positive S^r) has 
also been observed at even stronger condensates. 

We observe that the condensate formation substan- 
tially increases flatness: from about the Gaussian value 
of 3 for the weak condensate (solid squares in Fig. |2Fb)) 
to F ~ 5.5 for the strong condensate (solid squares in 
Fig. RTf)). This can be explained as follows: conden- 
sate shear suppresses turbulence level by stretching and 
destroying vortices. Note that strong vortices with the 
vorticity exceeding the external shear survive. Therefore, 
strong fluctuations are affected by condensate less than 
the mean level, which increases F. 

It is important to note that similarity of our spectra 
to those of [HI [T3] does not necessarily mean that fc -3 
spectrum at large scales in the Earth atmosphere is also 
fed by the inverse cascade. To establish whether this is 
the case, one needs to analyze the atmospheric data in 
the way described here: subtract the coherent flow, recal- 
culate the second and the third moment of fluctuations 
and use Eq. [T] It is likely that the baroclinic (large-scale) 
instabilities play a role in forcing the large-scale flows. To 
model an external large-scale forcing in our experiments 
we added a large magnet on top of the small-scale forcing 
(as described in [B]) and found that the modifications in 
S3 are similar to those when the large-scale flow is formed 
via spectral condensation. The mean subtraction recov- 
ers the energy flux from small to large scales in both 
cases. Similarly, the mesoscale turbulence in the Earth 
atmosphere should be affected by the large-scale flow re- 
gardless of its origin. Let us stress that our experimental 
system is much simpler than the Earth atmosphere. With 
regard to the meso-scale and large-scale atmospheric mo- 
tions, there are two most important differences, namely 
the character of stratification and the absence of rotation 
in our system. Recent numerical simulations (see [23] and 
references therein) show that stratification may enforce 
a 3D dynamics and the forward energy cascade. On the 
other hand, recent experimental studies of decaying tur- 
bulence suggest a strong role of rotation in establishing a 
quasi-2D regime in which geostrophic dynamics is dom- 



inant and the energy cascade is inverse (regime of low 
Froude and Rossby numbers) |24j . More experiments in 
forced turbulence are needed to understand the compet- 
ing effects of rotation and stratification along with the 
complex interplay between turbulence and waves, reso- 
nant wave- wave interactions, etc. True nature of atmo- 
spheric turbulence both at large and meso-scales can only 
be revealed by the atmospheric measurements. What we 
have shown here is the need to separate mean flows and 
fluctuations to recover the energy flux. 

We conclude by stating that the condensate strongly 
modifies both turbulence level and its statistics; different 
velocity moments are affected at different scales. 
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